Isospin-symmetry-breaking corrections to superallowed Fermi (3 decay: 
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We study the formalism to include isospin-symmetry-breaking corrections when extracting the up- 
down Cabibbo-Kobayashi-Maskawa matrix element from superallowed + — > + nuclear /3 decay. 
We show that there are no first order isospin-symmetry-breaking corrections to the relevant nuclear 
matrix elements. We find corrections to the treatment of Towner and Hardy, and assess these using 
schematic models of increasing complexity. 
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I. INTRODUCTION 

Superallowed Fermi (3 decay provides the most strin- 
gent test of the conserved-vector-current (CVC) hypoth- 
esis, the most precise value for the up-down Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element V u d, and the 
best limit on the presence of scalar interactions. With 
the confirmation of CVC, V u d can be extracted with 
great precision to test the Standard Model [H, 0, HI ■ For 
this, one needs to evaluate ~ 1% theoretical corrections 
that arise due to nucleus-dependent isospin-symmetry- 
breaking (ISB) effects between the parent and daughter 
states, and due to radiative effects 0, Q. These cor- 
rections are small, but significant, and their associated 
theoretical errors at present dominate the uncertainty of 
V u d because of the very high precision reached experi- 
mentally 0. 

In the 2005 survey of Hardy and Towner the re- 
sults for the set of superallowed + — > + transitions were 
statistically consistent, after including these theoretical 
corrections. However, Penning-trap measurements of the 
transition energy for 46 V @, [1] moved this case to more 
than two standard deviations away from the 2005 survey. 
This lead Towner and Hardy (TH) @ to reexamine their 
treatment of ISB corrections and to include the contribu- 
tion from core orbitals. The latter were found to be espe- 
cially important for 46 V and this anomaly disappeared. 

In this paper, we study the formalism to include ISB 
corrections, and contrast the TH treatment to exact re- 
sults. Before proceeding, we review the necessary theo- 
retical background, following the discussion in TH 

Superallowed + — > + Fermi (3 decay depends only on 
the vector part of weak interactions, and with CVC, the 
decay transition 11 ft value" should be independent of the 
nucleus: 
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where Gy is the vector coupling constant and Mf is the 
Fermi matrix clement. CVC depends on the assumption 
of isospin symmetry, which is not exact in nuclei, but 
broken by electromagnetic and quark mass effects. As a 
result, M F is reduced from its symmetry value of M = 
V2 for T = 1 parent and daughter states. Following TH, 
we introduce the ISB corrections 5c to the Fermi matrix 
element by 
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In addition, there are radiative corrections to Eq. |T]), 
but we focus on Sq here. These isospin corrections are ~ 
1%, but must be calculated with a theoretical uncertainty 
of 10%, to guarantee a desired accuracy of 0.1%. This 
presents a challenge for nuclear theory. 

Hardy and Towner have shown [l], Q that the calcu- 
lated corrections eliminate much of the considerable scat- 
ter present in the uncorrected ft values, and the sta- 
tistical consistency among the corrected values is evi- 
dence that the corrections have been reasonably com- 
puted. However, the importance of precisely testing the 
Standard Model stimulates us to undertake a reevalua- 
tion. With this, we wish to start and stimulate further ef- 
forts to systematically improve ISB corrections, based on 
an accurate understanding of ISB in nuclear forces [oL flp|| . 

This paper is organized as follows. In Sect. [Til we 
show that TH do not use the isospin operator to cal- 
culate 8c (as mandated by the Standard Model). To 
examine potential consequences of this, we review the 
TH treatment in Sect. IIIII A complete formalism is pre- 
sented in Sect. IIVI where we show that there are no first 
order ISB corrections to the relevant nuclear matrix ele- 
ments, which is also true for the work of TH. In Sect. W\ 
we compare the TH treatment to exact model results of 
increasing complexity, which can guide future improve- 
ments. We conclude in Sect. IVII 



II. 



TOWNER AND HARDY APPROACH TO 
ISB CORRECTIONS 
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In nuclei, the matrix elements of weak vector interac- 
tions are not modified by nuclear forces, except for cor- 
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rcctions due to ISB effects. Therefore, one has to eval- 
uate the contributions from electromagnetic and charge- 
dependent strong interactions to the Fermi matrix ele- 
ment Mp = (/|r+|i) between the initial and final states 
for superallowed (3 decay, \i) and |/), respectively. Here 
t+ is the isospin raising operator. 

Towner and Hardy [6( use a second quantization for- 
mulation to write the Fermi matrix clement as 

M F = Y,(f\a>U \i)(a\T + \l3) ) (3) 

where creates a neutron in state a and ap annihilates 
a proton in state (3. Thus, the label a is used to denote 
neutron creation and annihilation operators, while (3 is 
used for those of the proton. This notation is different 
from the standard notation [Til ] . in which b a is used to 
denote proton annihilation operators. 

The single-particle matrix element (a|r + |/3) is assumed 
to be given by the expression 

/>oo 

Hr+l/3) = 5 a ,f, K(r)R p Jr)r 2 dr = S a ,pr a , (4) 
Jo 

where R£{r) and R^(r) are the neutron and proton ra- 
dial wave functions, respectively. The problem is that the 
correct superallowed beta decay operator in the Standard 
Model is the plus component of the isospin operator. The 
operator in Eq. ([3]) is not the isospin operator, because 
the states \a) and \/3) are not the same. Instead, t + of 
Eq. is theplus component of the W-spin operator of 
MacDonald [ijj , which is reviewed in Ref . [ll[ . In addi- 
tion, Eq. © assumes that the radial quantum numbers 
of the states a and (3 must be the same. This need not 
be so. As a result, the Standard Model isospin commu- 
tation relations maintained in the W-spin formalism are 
lost. 

To obtain the commutation relations, we observe that 
Eqs. © and ^ correspond to the second-quantized 
isospin operators 

t + = E Sa <p Ta a ^ a p > ( 5 ) 

a, fj 

ct,j3 

so that 

[r+,r_] =^]M 2 ala a - ^ M 2 a^ap ^ r . (7) 

a (3 

The Standard Model isospin commutation relations are 
violated if one uses the isospin operators of TH. 

This formal problem motivates us to reevaluate the 
treatment of ISB corrections, and to study whether there 
are potential corrections to the extraction of V uc i- To this 
end, we review the details of the TH procedure for 5c- 
Although Eqs. ([3]) and ([J| are not formally correct, they 
do account for the important correction: the effects of 
the Coulomb interaction on the radial wave functions. 



III. TH TREATMENT OF 5c 

Towner and Hardy Q proceed by introducing into 
Eq. © a complete set of states for the (A — l)-particle 
system, \tt), which leads to 

M F = J2(f\^)(^K\i)rl. (8) 

a,7r 

The TH model thus allows for a dependence of the radial 
integrals on the intermediate state w. 

If isospin were an exact symmetry, the matrix elements 
of the creation and annihilation operators would be re- 
lated by hermiticity, (7r|a a |z) = (/|aj,|7r)*, and all radial 
integrals would be unity. Hence the symmetry-limit ma- 
trix element in this model is given by 

M =£|(/|4|7r}| 2 . (9) 

Towner and Hardy divide the contributions from ISB into 
two terms. First, the hermiticity of the matrix elements 
of a a and a' a will be broken, and second, the radial in- 
tegrals will differ from unity. Assuming both effects are 
small, TH calculate the resulting ISB corrections as @ 

5c = Sci + S C 2 , (10) 

where in evaluating 5ci ah radial integrals are set to unity 
but the matrix elements are not assumed to be related 
by hermiticity, and in evaluating 5c2 it is assumed that 
(n\a a \i) = (/|a|j7r)* but r£ ^ 1. We will study whether 
this is a useful representation of 5c- However, we em- 
phasize that the separation into 5ci and 5c2 is a model- 
dependent concept, inspired by the shell model Q. For 
example, this division is clearly model dependent when 
Mp is obtained from ab-initio calculations of the initial 
and final states, \i) and |/). In addition, we demonstrate 
below that this is not possible rigorously for schematic 
models. 



A. Radial overlap correction 5c2 

Towner and Hardy find that the radial correction, 5c2, 
is the larger of their two model corrections 0, H, Q . The 
Fermi matrix clement relevant for 5c2 is given by 

m f = EK/kk)! 2 ^, 

a,7r 




where f2J = (1 — r£) is a radial-mismatch factor. With 

the definition of the ISB correction factor in Eq. , TH 
approximate 5c2 by 

tf «« J-EK/I^)! 2 ^- (12) 
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Consequently, large contributions to Sc2 come with a 
large spectroscopic amplitude and a significant radial 
mismatch. 

In evaluating Sc2 of Eq. fT2")) . TH use guidance from 
experiment. Their results are based on shell model cal- 
culations of the spectroscopic amplitudes, but limit the 
sums over orbitals a and intermediate states it to those 
for which large spectroscopic factors have been observed 
in pick-up reactions. For 46 V, TH Q use this strategy to 
include two sd core orbitals, s-\_/2 and CZ3/2, in addition to 
the orbital of their earlier calculations. Their new 
result for &ci is 0.58% (see Table I in Ref. @), which is 
almost a factor two larger than the 2002 value Q . 

For the radial integrals, TH use the strong constraint 
that the asymptotic forms of all radial functions must 
match experimentally measured neutron and proton sep- 
aration energies. In many cases, TH have to truncate the 
model space to keep the calculations tractable. Their fi- 
nal values for Sci range between 0.17 and 1.50%, and 
increase with mass number (see Table II in Ref. [(|). 



two theorems, which show that there are no first-order 
ISB corrections to Fermi matrix elements. This provides 
a perturbative expansion, which allows for a simple esti- 
mate of 5c ■ 

We use the correct isospin operator 



(13) 



where a represents any single-particle basis, and cre- 
ates neutrons and b a annihilates protons in state a. The 
Fermi matrix element is then given by 



M F = (f\r+\i) 



(14) 



with |i) and |/) the exact initial and final eigenstates 
of the full Hamiltonian H = Hq + Vc, with energy Ei 
and Ef, respectively. Here Vc denotes the sum of all 
interactions that do not commute with the vector isospin 
operator T = J2t=i T i/2, 



B. Isospin-mixing correction Sci 



[H, T] = [V c , T] ^ and [H , T] = . 



(15) 



The isospin-mixing correction Sci is obtained by set- 
ting all radial integrals to unity, but including ISB cor- 
rections to the matrix elements of the creation and an- 
nihilation operators, (/|ajj7r)* and (7r|a Q |i). These arise 
because the neutron-rich and proton-rich states are dif- 
ferent. TH find that calculations of Sci are very sensitive 
to the details of the shell-model computation, but try to 
reduce the model dependence by using various experi- 
mental information [6j. 

To obtain Sci, TH use experimental single-particle en- 
ergies (on top of the core of the shell model calculation) , 
which differ for neutrons and protons. In addition, they 
include a two-body Coulomb interaction among the va- 
lence protons and increase all T = 1 neutron-proton ma- 
trix elements (relative to the neutron- neutron ones), so 
that the measured b and c coefficients of the isobaric 
multiplet mass equation (IMME) are reproduced. Fi- 
nally, TH account for weak transitions that can occur to 
non-analog + states. The adopted values for Sci range 
between 0.01 and 0.35%, and also increase with mass 
number (see Table III in Ref. @). 



IV. EXACT FORMALISM AND THEOREMS 
FOR ISB CORRECTIONS 

In this section, we present an exact formalism, inde- 
pendent of feasibility. We use this formalism to derive 



J 



We will use round bra and ket states to denote the 
eigenstates of the isospin-symmctric part of the Hamil- 
tonian, so Ho\n) = En\n). Obtaining the states \n) 
requires a solution of the ^4-body problem. 

The full initial cigenstate |z) can then be written as 



i) = VZi 



1 



MV C \i) 



Ei - KiHL, 
with projector A, = 1 — or equivalcntly, 

Similarly, the full final eigenstate |/) is given by 

1 



l/> = Vz. 



I/) 



E f -AfHA 



■A/Vc|/) 



(16) 



(17) 



(18) 



with A/ = 1— |/)(/|. The factors Zi and Zf are taken to 
be real and ensure that the full eigenstates \i) and |/) are 
normalized. As a result, it follows (due to the projection 
operators) that the deviations of Zi and Zf from unity 
start at second order in V c . 

We now evaluate Eq. (|T4"]) between the exact eigen- 
states given by Eqs. (JTOJ and (HI]). With (/|r+A 4 = 
and A/t + \i) = 0, we obtain 



Mp 



V Z i Z f 



M + (f\V c A f 



E f 



1 



AfHAf 



Ei - A,HA,, 



AiV c 



(19) 
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where M = (f\r+\i). Since Z itf = 1 + 0(V 2 ), it follows 
that ISB contributions start at second order. This is our 
first theorem and demonstrates that there are no first- 
order ISB corrections to Mp. 

We obtain a simpler form by expanding in the dif- 
ference of the charge-dependent interactions AVc be- 
tween the initial proton-rich and final neutron-rich states. 
Hence, AVc includes all charge-dependent interactions 
of the extra proton with the other nucleons in the initial 
state. In this case, we have 



|/> = |/) and 



Z\i) + 



1 



Ei — A{HoAi 



A,AV C \i) 



(20) 



where the first expression defines AVc an d -ffo includes 
the effects of Vc common to the initial and final states, 
for example the Coulomb interactions in the core. 

In this case, the final state is an eigenstate of Ho and 
obeys (/| r+Aj = 0. As a result, it follows that 



M F = VZM 



(21) 



This is our second theorem. As already shown, there are 
no first-order ISB corrections to Fermi matrix elements, 
and in this case 5c = 1 — Z has a straightforward per- 
turbativc expansion in AVc, starting at second order: 



1 



A;AV C \i)+2Re 



(i\ AVc As 



r 



i 



KAVc 



1 



Ei - A ? i7 A, 



AiAVfc \i) 



which follows from the normalization condition = 1. 
To second order in AVc, the full energy Ei can be taken 
as the ener gy e \ ^ of ffo- Examining the third-order term 
of Eq. (|22|) |17l | also shows that it is impossible to separate 
5c into two distinct terms. This is because it is not 
possible to distinguish whether the middle AVc is part of 
a correction to an intermediate state \ir) or to the initial 
state |z). Finally, we note that the two theorems are more 
general versions of the theorem of Behrends-Sirlin [l3[ for 
CVC in nucleons and of Ademollo-Gatto [l4| for weak 
decays of kaons. 



inert core of charge Ze, assuming harmonic-oscillator 
single-particle wave functions with oscillator frequency 
lu w 39McVA~ 1/3 . The nuclear Coulomb potential 
arises from the convolution of Ze 2 /(47r|r — r'Q with the 
charge density pc(r'). If we take the latter to be a con- 
stant within r R, the one-body Coulomb potential 
takes the simple form: 



AVc(r) 



Ze 2 
4ttR 



3 r 2 \ n , ,R 



(23) 



A. Simple estimate of Sc 

As an illustration of the above formalism, we cal- 
culate 5c for the case of a single particle outside an 



With R = 1.1 fm A 1 / 3 , we have Ruo = 0.22, independent 
of A, and therefore the correction scales as 5c ~ Z 2 . 
To make an estimate, we take the state \i) to be in the 
single-particle orbit with radial quantum number n = 
and angular momentum I. Using Eq. (|22[) . we find 



J 



^ ^ 1 1 — 



AUtt) 2 R 2 uj 2 ' n" 

v ; n>0 



r 2 drR i(r) 



r 2 \ 



etr -R)- 

r 



Rni{r) 



(24) 



r 



We calculate the summation numerically taking R equal 
to the oscillator length. This leads to 



5 C (0) = 0.0020% Z 2 
5 C {2) = 0.00071% Z 2 



5 C (1) = 0.0013% Z 2 
5 C (3) = 0.00043% Z 2 . 



(25) 



For a Z = 20 core and I = 3, we find <5 C (3) = 0.17% and 
thus S c = 3.33 x (5 C (3) = 0.57%, where the factor of 3.33 
arises from the Macfarlane-French sum rule UM for the 



three protons in the fj/2 orbit. This result is in qualita- 
tive agreement with the TH contribution of 0.45% (see 
Table I for 46 V in Ref. 0|). This indicates that Eq. |[2"T]) 
could be a useful starting point for realistic calculations. 
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V. SCHEMATIC MODELS 

Next we present exact evaluations of the Fermi matrix 
element for schematic models of increasing complexity, 
and compare our results to the treatment of TH. 



A. One-body problem 

Consider starting from the exact formalism. We can 
derive an effective single-particle potential U + Uc, where 
Uc accounts for charge-dependent effects and acts only 
on protons. The single-particle potential is introduced to 
minimize the effects of residual interactions 



and excited state. Then the exact eigenstate is a two 
component wave function, where the upper component 
represents the single particle plus unexcited core, and 
the lower component has the core in the excited state. 
The core excitation need not have angular momentum 
J = 0, but the coupled state is + with T = 1. In this 
case, the Hamiltonian is given as a two-by-two matrix: 



Hn = 



H sp AV 
AV H SD 



(32) 



where the second, lower component has a higher single- 
particle energy than the upper component. The eigen- 
functions are given by two-component "spinor" wave 
functions, for example for the initial state: 



AV = V + V c - (U + U c ) 
Then the Hamiltonian is given by 



H = T + U + U c + AV = H, 



sp 



AV. 



(26) 



(27) 



a-i Ui(v) 



(33) 



with normalizations given by 



In the simplest case, we assume that the one-body Hamil- / dr\Ui(r)\ 2 = / drlLJi 
touia.n ff™ is dominant. Thus we take the initial a.nH final J J 



= 1 and 



tonian H sp is dominant. Thus we take the initial and final 
states to consist of a single nucleon outside an inert core 
|0) . The core and nucleon have quantum numbers so that 
the coupled state is + with T = 1. 

The one-body basis states can be taken as eigenstates 
of the full single-particle Hamiltonian T + U + Uc, which 
we denote by |a), or by the eigenstates \a) of the isospin- 
symmetric part T + U. Here Uc is the difference between 
the proton and neutron potentials, which corresponds to 
Uc = Up — U n + Vc in the TH notation. The parent and 
daughter states are then given by 



PI 



(34) 

and where we have taken a, (3 to be real for simplicity. 

The presence of the charge-dependent interaction Uc 
(in H sp ) and of Vc (in AV^) causes the initial and final 
state values of a, (3 and their radial wave functions to dif- 
fer. In this model, the single-particle wave functions for 
i and / represent directly the single proton and neutron. 
The exact value of Mp is thus given by: 



|i) = 6t|0) and |/)=4|0) 



(28) 



It is convenient to express the isospin raising operator 
t + in a mixed representation. The creation operators of 
the two bases are related by 



M F /M = a f a z / dr U}(r) U^r) 

+ P f Pi [ drL* f (r)Li(r), (35) 



E4 



a\a) 



(29) 



since the core and its excitation are orthogonal. This 
may be rewritten as: 

Mf/Mq - i = -«? nto - 01 n< 2 > 



and therefore the isospin operator of Eq. p3j) reads 

r + = E4(«i«)^- ( 3 °) 

a, a 

This equation leads to an expression for Mp that is very 
similar to Eq. ([4]) of TH with the important difference 
that the states a and a need not have the same radial 
quantum numbers: 



+ (a f - ai ) a (i - nW) + (i3 f - ^ (i - n^) 



where in the TH notation 



= 1- / drU}(r)Ui(r), 
= I- J dvL){v)U{v). 



(36) 

(37) 
(38) 



M F /M = / r 2 drR*~{r)R a {r). 



(31) 



B. One-body problem with a single core excitation 

The simplest generalization of the previous problem 
is to allow the core to have two states, a ground state 



We next use the strategy of TH to evaluate this two- 
state core model. The states \w) consist of the ground 
state core and its excitation labeled by 1,2. For each of 
these, there is only one value for the single-particle index 
a. Therefore, the TH result for this model reads 

Mj H /M -l w -a 2 SlW-$QW 

+ (a f -ai)ai + (J3f- A) A- (39) 
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The contributions on the first line of Eq. (|39|) correspond 
to 5q2 and the terms on the second line to 6ci ■ 

In comparison with the exact result, we observe that 
TH neglect terms of order (a/ — ai)^ 1 ' . The relevant 
radial integrals arc of infinite order in Uc, so that setting 
them to unity in evaluating the second line of Eq. ([39]) 
may be significant relative to the required accuracy, in 
particular if the neutron and proton separation energies 
are very different. In addition, this schematic model indi- 
cates that the normalization conventions of Eq. (|34[) are 
just a choice, so that the separation into the two terms 
6ci, Sc2 seems rather arbitrary. It is just as reasonable 
to use the product aiUi(r) as the upper component Ui(r) 
and (3iLi(r) as the lower component Li(r) with the nor- 
malization J dr (|[/i(r)| 2 + |Z/.;(r)| 2 ). For that convention, 
the factors Q-ij, A./ would disappear from the formalism, 
and the separation of 5c into 5ci, Sc2 would neither be 
necessary nor possible. 



1. Evaluation using simple interactions 

Let the state |/) be governed by a harmonic oscillator 
Hamiltonian H , 



rr P 2 , 1 2 2 

Hn = — mu> r , 

U 2m 2 



(40) 



and the state |z) by the Hamiltonian H, with 

H = H Q + V C . (41) 

We use the inner form of the Coulomb potential, Eq. (|23p. 
to obtain simple expressions that show the order of vari- 
ous terms, 



Ze 2 r 2 
Vc(r) = -— — 
C( 1 4ttR2R 2 



(42) 



This is a qualitative approximation that has been tradi- 
tionally used to assess the size of various effects (Tlj . 

In this case, the Coulomb interaction shifts the square 
of the oscillator frequency from ui 2 to oj 2 (1 — 6), where 
the shift 6 is given by 



6 = 



Ze 2 



ATrmuj 2 R 3 



(43) 



We will consider the two lowest states with n = 0, 1 and 
angular momentum I = to study the effects of configu- 
ration mixing. 

Then for the final state single-particle basis, we have 



(r|0) 
(r|l) 



1 . e -r 2 /2b 2 



(7r&2)3/4 



3b 2 J (tt& 2 ) 3 / 4 



(44) 
(45) 



with oscillator length b — (muj) x l 2 . The initial state 
wave functions (r|0) and (r 1 1) have the same form, but 
with oscillator length 



b 2 



1-5/2' 

With this, we find the radial overlaps: 



(0|0) 
(1|1) 

(o|T) 
(i|o) 



1 - 6/2 



(1-6/4) 2 



3/4 



= 1 - 



3f_ 
64 



(46) 



0(6 3 ), (47) 



(1 - £/2) 3 / 4 (1 - 6/2 - 36 2 /32) 

(1- (5/4)7/2 

VgjCL -J/2) 3 / 4 
(2 -(5/2) 5 /2 ' 

-(0|T). 



(48) 



(49) 



(50) 



Now consider configuration mixing between the states 
and 1 by strong interactions. As as a result, we have 



I/) 



a\O)+0\l), 
5|0) + j9|T), 



(51) 
(52) 



where we take a, (3, a and /3 to be real for simplicity, and 
the quantities a — a and (3 — 13 are of order <5 2 . 
The Fermi matrix element is then given by 

M F /M = aa (0|0) + (3(3 (1|T) + [a(3 - (3a\ (0|T) . (53) 

The first two terms are equal to 1 + 0(<5 2 ). Since a = a + 
0{6 2 ) and similarly for (3, the leading first-order part of 
(0|1) (see Eq. ([15])) thus cancels exactly, and the last two 
terms of Eq. (|53")l start at order 6 3 . This validates that 
there are no first-order ISB corrections to Mp, and shows 
that certain approximations can violate our theorems. 



C. Two interacting nucleons outside an inert core 

Next we assume an inert core and two interacting nu- 
cleons with total angular momentum J = 0. The initial 
state wave function can have components spread over dif- 
ferent single-particle configurations with radial, orbital 
and total quantum numbers n, l,j: 

\i) = ^2 (j> m ,j,- m \°> Q ) 

m ,ri2 ,l,j,m 

>< AT' nM bl lljtm bl 2l ._ m \0) . (54) 

For clarity, we have taken two protons (fot) 2 on top of a 
+ core, coupled to J = 0, Mj = 0. A similar expres- 
sion for the final state involves the amplitudes A™ 1 '™ 2 ' '\ 
which differ due to the effects of the charge-dependent 
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interactions. In this way, two-nucleon correlations are in- 
corporated in a limiting case of the formalism of Sect . lIVI 
The exact expression for Mp is then given by 



M F /M = A J 

ni ,n2 



l-n 2 ,l,j * ^™i,™2,£,j 



x / r 2 drRf l *(r)Kl ll (r), (55) 

with radial wave functions Rij- For example, n\ could 
correspond to states in the shell-model valence space and 
n' x to a high-lying shell due to strong interactions. 
The TH approximation for Mp would be 



and the TH approximation would be 



model space r 



m™/m « y: 

ni ,ri2 



ni,n 2 ,l,j\2 



1 - Q 



(Af 



ni,n 2 ,l,j* _ j^ni,n 2 ,l,j *\ ^n 1 ,n 2 ,l,j 



\B 



"•i ,n2,l,j |2 



i - n£ niU) ) 



I j^n 1 ,n 2 ,l,j * _ j^n 1 ,n 2 ,l,j *\ j^ni,n 2 ,l,j 



VI. CONCLUSIONS 



>iW)l 



(58) 



model space p 

M™/M M ^ |^ni,na.W|a ^ _ n («i*i)) 

ni.n 2 ,l,j 



(A 



ni,Ti3,J,3* _ j^n 1 ,n 2 ,l,j *\ A n!,n 2 ,l,j 



f 



A 



(56) 



We compare the exact result of Eq. (|55|) with the TH 
approximation Eq. (|5C 



1. We find corrections to the radial overlaps, because 
the quantum numbers n\ and n[ need not be equal. 

2. The exact result mixes in higher- lying configura- 
tions that are not within the TH model space. To 
incorporate ISB effects due to higher-lying states, 
one needs to evaluate their contributions to charge- 
dependent effective interactions. In particular, an 
interesting topic for future study is the renormal- 
ization from long-range Coulomb effects. 

3. As in the previous models, the radial integrals of 
Eq. (|55|) are of infinite order in Vc , so that setting 
them to unity in evaluating bci might not be very 
accurate. 



D. Two nucleons with a single core excitation 

This model combines those of the two previous subsec- 
tions. The core can be excited so that the two interacting 
nucleons are outside a core in its ground or excited state. 
Using the previously adopted notation, the exact value 
of Mp reads 



We have studied the formalism to include ISB correc- 
tions to Fermi matrix elements, motivated by the recent 
experimental achievements on + — > + nuclear (3 decay 
and by the work of Towner and Hardy 0] . This is a key 
challenge for nuclear theory and pivotal for extracting the 
up-down CKM matrix element, which provides precision 
tests of the Standard Model. 

We have shown that TH do not use the isospin operator 
of the Standard Model to calculate ISB corrections. It is 
also true that their separation 6c = 8ci + <5c2 is model 
dependent Using a complete formalism, we derived 
two theorems, which demonstrate there are no first-order 
ISB corrections to Fermi matrix elements. 

Towner and Hardy correctly include the leading part 
of the Coulomb effects on the radial wave functions. We 
have found corrections to the TH treatment, and con- 
trasted these to exact results in schematic models of in- 
creasing complexity. One of the differences is that the 
radial overlaps need not have the same radial quantum 
numbers (as assumed in TH). This mixing has also been 
pointed out in density-functional based calculations of 
ISB corrections In addition, significant to the re- 

quired accuracy, there may be contributions from higher- 
lying configurations that arc outside the model space. 
This requires a careful inclusion into charge-dependent 
effective interactions, where our accurate understanding 
of isospin-symmctry breaking in nuclear forces can be 
very helpful, and a careful study of truncation effects, 
where modern methods can lead to improvements. 

Numerical evaluations using the formalism presented 
here are needed as a next step. We hope that our work 
stimulates further efforts to make systematic improve- 
ments to the important problem of ISB corrections to 
superallowed transitions. 



M F /M = 

Til ,71^ ,7X2 



A n u n 2 ,l,j* A n lt n 2 ,l,j / ^ * (r) 



. (57) 
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